Heavy quasiparticles in the ferromagnetic superconductor ZrZn 2 
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We report a study of the de Haas-van Alphen effect in the normal state of the ferromagnetic su- 
perconductor ZrZn2. Our results are generally consistent with an LMTO band structure calculation 
which predicts four exchange-split Fermi surface sheets. Quasiparticle effective masses are enhanced 
by a factor of about 4.9 implying a strong coupling to magnetic excitations or phonons. Our mea- 
surements provide insight in to the mechanism for superconductivity and unusual thermodynamic 
properties of ZrZn2. 
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The recent reports of the coexistence of ferromag- 
netism and superconductivity in UGe2 Q, ZrZn2 0, 
URhGc j| as well as the discovery of superconductiv- 
ity in the non-magnetic hep phase of Fe under pressure 
Pj have reopened the debate regarding the relationship 
of magnetism and superconductivity. Motivated by the 
case of liquid 3 He, the search for superconductivity me- 
diated by spin fluctuations in nearly magnetic materials 
has a long history [fjj. The observation of superconduc- 
tivity in the ferromagnetic phase of ZrZn2 Q (as well as 
UGe2 and URhGe) together with the extreme sensitivity 
of its occurrence to sample purity |2[ |(| , make it a strong 
candidate for magnetically mediated and spin-triplet su- 
perconductivity, as is believed to exist in Si'2Ru04 || 
although other scenarios are not excluded . 

A detailed knowledge of the electronic structure of 
ZrZii2 is crucial to the understanding of both its normal 
and superconducting properties. Thus we report here 
a detailed study of the Fermi surface (FS) made using 
angle-resolved measurements of the de Haas- van Alphen 
effect (dHvA). Our results are compared with an ab-initio 
electronic structure calculation. We observe much of cal- 
culated Fermi surface which consists of four exchange- 
split sheets. Our study shows that ZrZn2 is character- 
ized by a large quasiparticle density-of-states (DOS) at 
the Fermi energy in the ferromagnetic state, which arises 
partly from the band structure and partly from a large 
mass enhancement. 

ZrZn2 has many unusual properties. Below Tfm = 
28.5 K it becomes ferromagnetic with an ordered mo- 
ment of 0.1 1[Ib per formula unit at low temperatures 
(T = 2 K). The intrinsic moment is unsaturated, with an 
applied magnetic field of 6 Tesla causing a 50% increase. 
Compared with other d-band metals, it has an extremely 
large electronic heat capacity at low temperatures C /T 
= 47 mJK~ 2 mol -1 . ZrZn2 crystallizes in the C15 cubic 
Laves structure, with lattice constant a — 7.393 A, the 
Zr atoms forming a tetrahedrally co-ordinated diamond 
structure. Its magnetic properties derive from the Zr Ad 
orbitals, which have a significant direct overlap. 

Interest in ZrZn2 has been rekindled by the discovery 
that it is superconducting at ambient pressure and that 
both the superconductivity and ferromagnetism are si- 



multaneously destroyed by the application of hydrostatic 
pressure. For the samples considered here, the onset 
temperature for the superconductivity is Tsc ~ 0.6 K 
with B C 2 ~ 0.9 T. The most striking feature of ZrZii2, 
UGe2 and URhGe is perhaps that the same electrons 
are thought to participate both to superconductivity 
and ferromagnetism, in contrast with the situation in 
other "magnetic" superconductors e.g. borocarbides, 
RuSr2GdCu208 @], where the magnetism and supercon- 
ductivity occur in different parts of the unit cell. Further- 
more, ZrZn2 is a three-dimensional intermetallic com- 
pound in which the itineracy of the d-electrons is almost 
unquestionable, whereas some doubts remain about that 
of the 5/ electrons and roles of the strong magnetocrys- 
talline anisotropy and quasi-two-dimensional electronic 
structure in UGe2 and URhGe. 

The dHvA effect is due to the quantization of the cy- 
clotron motion of charge carriers and results in a mag- 
netization, M osc , which oscillates with magnetic field, B. 
We interpret our results using the Lifshitz-Kosevich ex- 
pression [[[o| in which each extremal area A of a FS sheet 
perpendicular to B gives a contribution to the oscillatory 
magnetization, 
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where F(B) = (h/2ne)A(B). R D and R T are the Dingle 
and temperature damping factors, respectively. The tem- 
perature dependence of M osc yields the cyclotron mass, 
to*, through R T = Xj sinhA where X — 27r 2 fc B T/7iw c 
and uj c = eB/m*. This mass includes all many-body 
renormalizations and is related to the corresponding or- 
bit through, 
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27T JFS orbit VF 

The quasiparticle lifetime, r, and the mean-free-path, I, 
can be determined on each orbit via the ficld-dcpcndcnt 
factor Rd = exp(— irmb/eBr). Eq. (|l]) shows no "spin- 
splitting factor" [|l0| since in a ferromagnet, spin-split 
sheets of the FS are resolved separately. 

The dHvA magnetization was measured using a stan- 
dard low-frequency field-modulation technique. Data 
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were collected using a top-loading dilution refrigerator 
equipped with a 13.5 T magnet reaching temperatures 
as low as 30 mK and with a 3 He system equipped with 
an 18 T magnet. High-quality samples (p =0.6 fiflcm) 
were prepared by directional cooling of a ZrZn2 melt con- 
tained in an Y2O3 crucible inside a tantalum bomb [ pT| . 

We performed angle-resolved dHvA measurements for 
fields applied perpendicular to the [001] and [110] di- 
rections. Fig. 0(a) shows a typical Fourier transform of 
our dHvA data collected for B close to [111]. The spec- 
trum is dominated by the low frequencies F a — 1.42 kT 
and Fp — 1.60 kT and their harmonics. These frequen- 
cies have been previously assigned to an exchange-split 
pocket of the Fermi surface centered on the V point of the 
Brillouin zone [see Fig. 0(a)]. The frequency difference 
F a — Fp is a measure of the ferromagnetic exchange split- 
ting in ZrZn2. In addition to a and (3, we also observe 
a number of closely spaced frequencies in the range 1.6 
- 2.5 kT (77, k, 9, A, /i, and 7) and a strong component 
(F v ) near near 3.2 kT. Some frequency branches (/i and 
7) appear as multiplets. In all, our study revealed 17 
frequency branches, whose angular dependence is shown 
in Fig. |l|(b) . Our results are in good agreement with pre- 
vious dHvA studies of ZrZn2 in which the a, (3 and v 
branches were observed [Q . In addition to extremal FS 
areas, the temperature dependence of the dHvA ampli- 
tude enables the determination of the cyclotron effective 
mass corresponding to each orbit (see Table |) from the 
factor R T in Eq. (§). 

In order to understand our results we performed self- 
consistent spin-polarized band structure calculations us- 
ing the LMTO method within the local-spin-density ap- 
proximation (LSDA) We obtained a relaxed lat- 
tice constant, a — 7.223 A and a magnetic moment, 
(i = 0.18 ^ts/Zr in excellent agreement with the mea- 
sured one (of 0.17 yUs/Zr Q). Our calculation agrees 
well with others in the literature |8|, [l4| . The resulting FS 
and DOS are shown in Fig. ^| The angular dependence 
of the dHvA frequencies [see Fig. |l](c)] imposes strong 
constraints on the possible topologies of the FS and is 
therefore an important tool in verifying the validity of 
the calculated FS. The cyclotron band mass for each ex- 
tremal orbit (see Table ||) is obtained from the relation 
rrib = (h 2 /2tx) (dA/dE) Ep by numerical differentiation 
with a AE of 0.5 mRy 

The ab-initio calculated FS appears to reproduce the 
essential topology of the observed one. We assign a and 
P to the spin-split spherical FS sheet as in previous 
studies The complex of branches observed near [111] 
including 7, 77, 9, k, A and p appears to be due to orbits 
associated with the X-centered "pillow" features of band 
27. Both spin \ and | FS sheets show this feature, thus 
there are pairs of exchange-split branches such as 9 and 
7 and (rj, n) and (A, /1). The 6*-branch is probably ob- 
scured by the strong a branch near the [110] direction, 
only allowing detection near [111]. Another prominent 



TABLE I: Experimental and calculated dHvA frequencies 
and cyclotron masses for the principal high-symmetry direc- 
tions in ZrZn2 for (B) — 9.2T. Branch assignments refer to 
Fig. jij(b) and the FS orbits are denoted by orbit center, band 
index and spin (| is the majority). 





Experiment 




Calculation 






branch 


F (kT) 


m* / m e 


FS Orbit 


F (kT) 


mi / m e 


[110] 


CX 


1.425 


0.81(4) 


^30, J. 


(sphere) 


1.193 


0.24 





1.600 


0.95(5) 


1^30, T 


(sphere) 


1.329 


0.26 


7 


1.810 


1.27(3) 


X 27,i 


(pillow) 


1.438 


-0.51 


6 


3.400 


3.7(10) 


X.28.T 


(dog bone) 


4.227 


-2.19 


€ 


4.770 


3.9(2) 




(dog bone) 


5.538 


-1.53 


c 


5.300 


3.5(10) 
















X-29.T 


(dog bone) 


5.541 


-4.25 


[111] 


a 


1.425 


0.82(5) 


r 30,i 


(sphere) 


1.193 


0.24 





1.600 


1.00(4) 


1^30 ,| 


(sphere) 


1.328 


0.26 




1.740 


3(1) 




(pillow) 


1.432 


-1.16 





1.809 


2(1) 


X 2 7,T 


(pillow) 


1.376 


-0.64 


K 


1.860 


3(2) 


X27.T 


(pillow) 


1.376 


-0.64 


Ai 


2.140 


1.6(1) 


X 27.i 


(pillow) 


1.833 


-0.76 


A 2 


2.180 




X 27,i 


(pillow) 






7i 


2.240 


1.9(1) 


X 27,i 


(pillow) 


1.833 


-0.76 


72 


2.272 


3.1(2) 


X 27,i 


(pillow) 






Mi 


2.307 


1.4(3) 


X 27,i 


(pillow) 


1.833 


-0.76 


A12 


2.330 


2.4(4) 


X 27,i 


(pillow) 






"i 


3.170 


2.20(6) 




(neck) 








3.230 


1.60(3) 


L 28,i 


(neck) 


3.067 


0.53 


"3 


3.260 


1.45(5) 




(neck) 












L29,T 


(neck) 


2.250 


7.89 








L28, T 


(neck) 


3.517 


1.33 


[001] 


? 


0.645 


4(2) 


W 2 7.i 


'pillow-neck) 


0.767 


0.68 


a 


1.420 


0.9(1) 


r 30,i 


(sphere) 


1.201 


0.24 





1.610 


1.1(1) 


r30 ? f 


(sphere) 


1.340 


0.26 


7T 


1.970 


4(1) 


W28,T 


(rosette) 


2.626 


-1.28 


P 


2.270 


7(2) 


W 2 8.i 


(rosette) 


3.232 


-0.94 


a 


3.630 


5(2) 


W 29 , T 


(rosette) 


3.699 


-7.96 



feature of the rotation plot is the v branch which has a 
minimum frequency for B || [111]. Some branches, e.g. \i 
and 7 [Fig. 0(a)], show small splittings which could arise 
from extremely small corrugations in the FS or a small 
misalignment of the sample. The v branch is a signature 
of the neck around the L point in the [111] direction and 
is in excellent quantitative agreement with the neck for 
the minority spin (|) of band 28. Although three sheets 
of the FS are predicted to have necks along the [111] 
directions (orbits L28,j,, L28,|, 1^29^) only one neck-like 
branch is observed experimentally. The most likely rea- 
son is that the higher masses of the L28,! anc ^ -^29,1 orbits 
(see Table |) make the signals too weak to be resolved. 

The topology of sheet 28 (J, and \ spin) is reminiscent 
of the FS of Cu ]l5] ]. Thus, we expect to see the other 
orbits seen in Cu. The S and e orbits near B || [110] 
(X 2 8,| and X 2 8,J correspond to the "dog bone" in Cu. 
We identify the ir and p branches as the W28,i an( ^ Was.T 
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FIG. 1: (a) A dHvA spectrum collected for B applied 1.5 deg from [111] towards [001]. Data was collected in the range 
7 < B < 13.4 T at T=25 mK. (b) The angular dependence the dHvA frequencies, (c) Calculated dHvA frequencies. 



calculated orbits and as the "four-cornered rosette" or- 
bits in Cu for fields B || [001]Jl||. These orbits connect 
the [111] necks in neighboring Brillouin zones and are 
shown in Fig. |^(a). The low- frequency £ branch appears 
to be due to the Wn\- Finally, the a branch appears to 
be due to either the rosette-type orbit VF29,t of sheet 29 or 
an extension of the branch to [001] which may oc- 

cur with a small energy shift of band 27. Our calculation 
overestimates the size of the "dog bone" and "rosette" or- 
bits on sheets 28 and 29. We note that because all these 
orbits are related to the necks, both frequencies would 
be decreased if these bands were lowered relative to Ep 
or if the necks opened out more rapidly along [111] than 
the calculation predicts. 

We now consider the mass enhancement of the quasi- 
particles in ZrZn2. Comparison of the measured and cal- 
culated cyclotron masses shows that the mass enhance- 
ments m* /nib = 1 + A are in the range of about 2-5 for 
B = 9.2 T. The linear term of the specific heat 7 shows 
a strong held dependence, 7 = 35 mJ mol _1 K~ 2 for B = 
9.2 T and 47 mJ mol^R- 2 for B = |@, thus we ex- 
pect m* to be field dependent. Comparing the measured 
7 for 5=9.2 T with 7ca i c = (tt 2 /3)k 2 B N tot (E F ) = 9.6 
mJ mol -1 K~ 2 |nj we find an average mass enhancement 
m* /mi, — 3.7 which is in the center of the range observed 
by dHvA. If we make the same calculation for B — 0, we 
find m'/nif, — 4.9. To our knowledge, no other alloy 
of 3ei-band metals shows a mass enhancements close to 
these. Palladium which is is often taken as an example of 
a system with strong spin fluctuations has an enhance- 
ment of about two. 

Table summarizes the important parameters for 



TABLE II: Selected measured (calculated) Fermi surface pa- 
rameters for ZrZn2 J20|. 



branch(orbit) 


a(r 30 ,i) 


/3(r 30 , T ) 


7(X 27 ,|) 


e(K 2 8,t) 


F (kT) 


1.42(1.19) 


1.60(1.33) 


1.81(1.38) 


4.77(5.54) 


(k F± ) (A- 1 ) 


0.21(0.19) 


0.22(0.20) 


0.23(0.20) 


0.38(0.41) 


m*,mb (m e ) 


0.82(0.24) 


0.95(0.26) 


1.27(-0.51) 


3.9(-1.53) 


up (Mms -1 ) 


0.30(0.92) 


0.27(0.90) 


0.21(0.45) 


0.11(0.31) 


r(ps) 


0.28 


0.31 


0.32 




l (A) 


2580 


2750 


1500 





some FS orbits in ZrZn2. The measured Fermi velocity 
and cyclotron mass include the effects of electron corre- 
lations. The quasiparticle lifetime and mean- free-path, I, 
are determined from the Dingle factor Rd- It is interest- 
ing to note that I ~ 1500-2800 A is considerably larger 
that the superconducting coherence length £ w 270 A in 
this superconducting sample allowing the possibility of 
unconventional superconductive pairing ||. 

It is interesting to compare the properties of ZrZn2 
with other d-band ferromagnets. Fig. ||(b) shows the 
paramagnetic DOS for ZrZii2. The Fermi level, Ep, lies 
in a region where the DOS is large due to the Ad Zr 
band. Also shown on Fig. ^(b) are the relative shifts of 
Ep for both spins in the ferromagnetic state. The small 
exchange splitting (A ks 70 meV) of ZrZn2 means that 
the Fermi levels for both spins remain in the high DOS 
region giving ZrZii2 its unusually large DOS in the fer- 
romagnetic state in contrast with Fe and Ni (Table |l|) . 
In the "stronger" ferromagnets, Nt ot (ep) is considerably 
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(a) 




E(eV) 



FIG. 2: (a) Calculated Fermi surface of ZrZri2 for the four 
bands crossing the Fermi level. The left and right parts of 
each figure are the majority (f ) and minority (J.) spin sheet 
respectively. The labeled solid lines are some of the orbits dis- 
cussed in the text, (b) The calculated DOS for paramagnetic 
ZrZn2 (see text). The vertical solid lines represent the energy 
shift of the spin-| and j Fermi energies in the ferromagnetic 
state. 



TABLE III: Comparison of band parameters (magnetic mo- 
ment, fi, exchange-splitting, A, DOS at Ef between ZrZn2 
and some d-band ferromagnets [n8| 





Fe 


Ni 


ZrZn 2 


H = n T - raj. {jis/tM.) 


2.12 


0.56 


0.16 


A (eV) 


2.2 


0.6 


0.07 


iV T (£ , F)(eV" 1 f.u.- 1 spin -1 ) 


0.83 


0.17 


2.39 


iVj( J BF)(eV- 1 f.u.- 1 spin" 1 ) 


0.25 


1.50 


1.80 


iVtot(^F)(eV- 1 f.u.- 1 ) 


1.08 


1.67 


4.19 



less in part due to the larger exchange splitting moving 
Ep for one of the spin out of the d-band. A high density 
DOS at the Fermi energy generally favors superconduc- 
tivity so this may be a key factor in understanding the 
superconductivity in this material. 

In summary, we have used angle-resolved dHvA mea- 
surements to characterize the Fermi surface of the ferro- 
magnetic superconductor ZrZn2. Our observations are in 
good agreement with an ab-initio LSDA band structure 
calculation and shed light on some of the unique proper- 
ties of ZrZii2. We find that ZrZn2 has a very large DOS at 
Ef compared to other d-band ferromagnets both because 
of the presence of several large sheets of FS and because 
the Fermi levels for both majority and minority spins lie 
in the Zr d-band in the ferromagnetic state. Further, 
the fermion quasiparticles are strongly enhanced with an 
average mass enhancement of about 4.9. This value is 
large for materials containing only d-band metallic ele- 
ments and comparable with that found in the oxide sys- 
tem Sr2Ru04 Jl^]. The presence of the high density of 
quasiparticles naturally explain the very large electronic 
heat capacity observed in ZrZn2 and also explains the 
unsaturated magnetic moment (or large intrinsic mag- 
netic susceptibility) in the ferromagnetic state. Finally, 
the high quasiparticlc DOS will favor superconductivity. 
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